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Abstract 



Parity and time invariance violating (P, T-odd) nuclear forces produce 
P, T-odd nuclear moments. In turn, these moments can induce electric dipole 
moments (EDMs) in atoms through the mixing of electron wavefunctions of 
opposite parity. The nuclear EDM is screened by atomic electrons. The EDM 
of an atom with closed electron subshells is induced by the nuclear Schiff mo- 
ment. Previously the interaction with the Schiff moment has been defined for 
a point-like nucleus. No problems arise with the calculation of the electron 
matrix element of this interaction as long as the electrons are considered to 
be non-relativistic. However, a more realistic model obviously involves a nu- 
cleus of finite-size and relativistic electrons. In this work we have calculated 
the finite nuclear-size and relativistic corrections to the Schiff moment. The 
relativistic corrections originate from the electron wavefunctions and are in- 
corporated into a "nuclear" moment, which we term the local dipole moment. 
For 199 Hg these corrections amount to ~ 25%. We have found that the nat- 
ural generalization of the electrostatic potential of the Schiff moment for a 
finite-size nucleus corresponds to an electric field distribution which, inside 
the nucleus, is well approximated as constant and directed along the nuclear 
spin, and outside the nucleus is zero. Also in this work the 239 Pu atomic EDM 
is estimated. 



PACS: 32.80.Ys,21.10.Ky,24.80.+y 
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I. INTRODUCTION 



The best limit on parity and time invariance violating (P, T-odd) nucleon-nucleon interac- 
tions (as well as quark-quark P, T-odd interactions) has been obtained from the measurement 
of the 199 Hg electric dipole moment (EDM) IJ. The mechanism of this EDM generation 
is the following. P, T-odd nuclear forces create P, T-odd nuclear moments, e.g. the EDM 
and Schiff moment (SM). According to the Schiff theorem the EDM of a point-like 

nucleus is completely screened by atomic electrons, so it cannot be measured. However, the 
electrostatic interaction between atomic electrons and the nuclear Schiff moment induces an 
atomic EDM. 

The electrostatic potential produced by the Schiff moment is usually presented in the 
form 

ip(R) =4ttS- V5(R), (1) 

where S is the Schiff moment (vector), <5(R) is a delta-function. The contact interaction 
— eif mixes s- and p-wave electron orbitals and produces EDMs in atoms; for example the 
atomic EDM induced in an atom with a single electron in state ns has the form 

^ (ns\ — eip{H)\mp) (mp\ — eH\ns) 

datom = 2 2^ —p ■ (2) 

The expression ([I]) is consistently defined for non-relativistic electrons. Using integration by 
parts, it is seen that the matrix element (s\ — ecp\p) is finite, 

(s\ — e(p\p) = AireS ■ (V^]^ p )r = o = constant . (3) 

However, atomic electrons near the nucleus are ultra-relativistic, the ratio of the kinetic 
or potential energy to mc 2 in heavy atoms is about 100. For the solution of the Dirac 
equation, (Vipt4 , P )R->o — ► oo for a point-like nucleus. Usually this problem is solved by 
a cut-off of the electron wavefunctions at the nuclear surface. However, even inside the 
nucleus Viplip p varies significantly, « Z 2 a 2 , where a is the fine-structure constant, Z is the 
nuclear charge. In Hg [Z = 80), Z 2 a 2 = 0.34. Recently, proposals have been made to 
measure EDMs of very heavy atoms like Ra [Z = 88) and Pu [Z = 94) jjlTJ where 

P, T- violating nuclear moments and the resulting atomic EDMs are very strongly enhanced. 

A consistent treatment of the Schiff moment is needed especially because the Schiff 
moment itself is defined as the difference of two approximately equal terms (see Eq. (|I5D). 
The aim of the present work is to develop a consistent theory of the nuclear Schiff moment, 
and the atomic EDM it induces, which properly takes into account the relativistic character 
of the electron wavefunctions inside the nucleus. 

For relativistic electrons we should introduce a finite-size Schiff moment potential. In 
this paper we show that the natural generalization of the Schiff moment potential for a 
finite-size nucleus is 

<p(R) = -^j^n(R) , (4) 

where B = J n(R)R 4 dR ~ R%/5, Rn is the nuclear radius, and n(R) is a smooth function 
which is 1 for R < R^ — 5 and for R > + ^ n(R) can be taken as proportional to the 
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nuclear density. This potential @ corresponds to a constant electric field inside the nucleus 
(see Fig. |X[) which can be produced by P, T-odd nuclear forces or by an intrinsic EDM of 
an external nucleon. This expression has no singularities and may be used in relativistic 
atomic calculations. 

A more accurate treatment requires the calculation of a new nuclear characteristic which 
we call the local dipole moment (LDM). This moment takes into account relativistic correc- 
tions to the nuclear Schiff moment which originate from the electron wavef unctions. So in 
the non-relativistic limit, Za -> 0, the LDM L — S. For 199 Hg, L « S(l-0.8Z 2 a 2 ) » 0.755. 
When considering the interaction of atomic electrons with the LDM we define it as placed 
at the center of the nucleus, that is the electrostatic potential is 

y?(R) = 4ttL • V5(R) . (5) 

This paper is organized in the following manner. In Section we derive a general 
expression for the dipole component of the P, T-odd electrostatic potential inside the nucleus. 



In Section |TTT] we take the electronic matrix element of this potential and show that it is 
related to the nuclear Schiff moment. In this section the electronic and nuclear problems are 
separated. It is shown that it is convenient to define a "nuclear" moment (the local dipole 
moment) which is the nuclear Schiff moment with higher-order (relativistic) corrections 
which originate from the electron wavefunctions. In Section [TV] we calculate various nuclear 
LDMs which arise due to P, T-odd nuclear forces; we calculate the contribution of an external 
proton and that of core protons to the LDM of a spherical nucleus, and we calculate the 
collective LDM of an octupole-deformed nucleus. Then in Section |V] we calculate the electric 
field distribution associated with the nuclear Schiff moment. In Section [VI] we estimate the 
size of the atomic EDM induced in 239 Pu. 



II. THE SCHIFF CONTRIBUTION TO THE NUCLEAR ELECTROSTATIC 

POTENTIAL 

The nuclear electrostatic potential with electron screening taken into account can be 
presented in the following form (see, e.g., Appendix of Ref. M for derivation): 

where ep is the nuclear charge density, / pd 3 r = Z, and d = e / prd 3 r = e(r) is the nuclear 
EDM. The second term cancels the dipole long-range electric field in the multipole expansion 
of y?(R). The Coulomb potential i^_ r t can be expanded in terms of Legendre polynomials 



^Pz(cos^), (7) 



|R-r| Y r > +1 

where r< (r>) is min[r, R] (max[r, R]). The P, T-odd part of the potential ([?]) originates 
from the odd harmonics I. The third harmonic I = 3 corresponds to the octupole field which 
has been considered in [11 1. The contribution of the I = 3 term is usually small (in 199 Hg, 



which has nuclear spin I — 1/2, it vanishes). Higher I always gives negligible contributions. 
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Therefore, we concentrate on I — 1 (it may be presented as ^-0(i? — r) + ^0(r — i?), 
where Q(r — R) = 1 for r > R and 0(r — i?) = for r < i?): 



rp(r)<i 3 r 



+ eR 



OO j. 

— p(r)d 3 r 
r r 6 



e(r) ■ R /•« 
Zi? 3 



p(r)d 3 r . 



(8) 



Note that in the second (screening) term in (§) we only keep the zero multipole / = 0. Also 
note that for R — > oo the first and third terms of Eq. (^j) cancel each other. Therefore, we 
can use f R = / °° — = — and present <^^^(R) as 



(1) (R) = eR- 



R 



Zi? 3 i? 3 



p{v)d 3 r 



(9) 



We see that ^^(R) = for R > R N (nuclear radius) since p(r) = in this area. We will 
see in the next section that this potential (|S|) is related to the Schiff moment. 

III. ELECTRON MATRIX ELEMENTS OF THE P, T-ODD ELECTROSTATIC 

POTENTIAL. 

All the electron orbitals for I > 1 are extremely small inside the nucleus. Therefore, we 
can limit our consideration to the matrix elements between s and p Dirac orbitals. We will 
use the following notations for the electron wavefunctions: 



V>(R) 



-i(<7 ■ n)g(R)Vtji m 



(10) 



where fijj m is a spherical spinor, n = H/R, and f(R) and g(R) are radial functions (see, e.g., 
fl~2|l). Using (cr • n) 2 = 1, then we can write the electron transition density as 



Psp {R) = ^ p = ntn p u sp (R) 

U sp (R) = f 8 (R)f p (R) + 9s(R)g P (R) = E£i hR k 



(11) 
(12) 



The expansion coefficients are calculated analytically and are presented in Appendix [A|; 
as is seen here, the summation is carried over the odd powers of k. Now we can find the 
matrix element of the electron-nucleus interaction, 



(R)b) = - e 2 ( S |nb)-{/ o °° [(i(r) 



r r 

U sp dR + 



, . U sp R dR 
r 6 Jo 



-e 2 (s 



n b) • { E 



:(r)(r fc+1 ) 



k + A 



(rr fe+1 )]} 



p d 3 r| 
(13) 



where (s|n|p) = / QlnQ p d<j) sin 6 d9 and (r n ) = J p(r)r n d 3 r. Note that all vector values (rr n ) 
are due to the P, T-odd correction to the nuclear charge density p, while j?{r n ) are the usual 
P, T-even moments of the charge density starting from the mean-square radius j?(r 2 ) = r 2 
for k — 1. 

In the non-relativistic case (Za — > 0) we have just b\ ^ 0, and so 
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e^\p) 



e 2 h 



(s\n\p) ■ — (r)(r 2 ) (rr 2 ) 



4vreS • (V^ P ) fl - 



•0 ) 



where the Schiff moment S is defined as 



e 
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(r 2 r) 



3Z 



(r 2 )(r) 



S I/I 



(14) 



(15) 



I is the nuclear spin. The expressions (14), (ll5f) agree with the results of Ref. 0. 

Therefore, Eq. (|T3p gives us the possibility of a consistent relativistic treatment of the 
atomic effects produced by P, T-odd nuclear forces. The nuclear and electronic problems 
can be separated in the following way. The nuclear calculations can provide us with the 
value of the local dipole moment (LDM) 



;& 1 (A; + l)(A; + 4) _ 



(rr k+1 ) 



k + 4 
3Z 



(r)(r fe+1 ) 



L I/I 



(16) 



which coincides with the Schiff moment S ( |i~5"|) in the non-relativistic limit (Za — > 0). Note 
that this "nuclear" moment contains relativistic corrections which arise from the electron 
wavefunctions which are calculated analytically. (It should further be noted that the cor- 
rections originating from the s-p\ji and s-p^/2 matrix elements are different.) The electron 
matrix elements are then given by 



(s\ - e<p w \p) = 4vreL • (V^ s V P )^o = 3eL • (s|n|p) 



' fsfp + 9s9p 
R 



(17) 



These formulae C |l6l , [TT|) , in principle, solve the problem of the consistent approach for the 
calculation of the interaction of the relativistic electrons with the Schiff moment. Note 
that to achieve ~ 10% accuracy it is enough to keep in L just the first correction, 63/61 = 
—3/5 Z 2 a 2 /R 2 N for the s-pin matrix element and 63/61 = —9/20 Z 2 a 2 /R 2 N for the s-p^/2 
matrix element (see Appendix [A]); and at this level of accuracy (~ 10%) the values of the 
coefficients 63/61 (for s-pi/2 and s-p^/2) can be taken to be the same. 



IV. LOCAL DIPOLE MOMENTS INDUCED BY P, T-ODD NUCLEAR FORCES 



We can now calculate local dipole moments induced by P, T-odd nuclear forces. We will 
begin in Section [IV A| with the calculation of the contribution of an external proton to the 
local dipole moment of a spherical nucleus. Because the best limit on the P, T-odd nucleon- 
nucleon interaction has been extracted from 199 Hg (which has an external neutron, so only 
the core protons contribute to the LDM) the result of Section [IV A| is not so interesting 
by itself. However, as will be explained in Section [IV B|, it provides us with a check of the 



method for the calculation of the contribution of core protons to the LDM. Then in Section 



IV Q we calculate the collective LDM of an octupole deformed nucleus. 



The P, T-odd nucleon-nucleon interaction, to first-order in the velocities p/ra, can be 
presented as || 
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W ab = 7/| ((VabVa - VbaO-b) ■ V a <5(r a - r 6 ) + r]' ab [cr Q x cr 6 ] • {(p a - p 6 ), 5(r a - r 6 )}) , 

(18) 

where { , } is an anticommutator, G is the Fermi constant of the weak interaction, m is 
the nucleon mass, and cr, r, and p are the spins, coordinates, and momenta of the nucleons 
a and b. The dimensionless constants r] a b and r]' ab characterize the strength of the P, T-odd 
nuclear potential (experiments on EDMs are aimed to measure these constants). 



A. Nuclear LDM produced by external proton 

In this section we calculate the LDM arising due to an external proton. We are therefore 
interested in the P, T-odd interaction of the external proton with the core nucleons. We 
can average the two-particle interaction fll8|) over the core nucleons to obtain the effective 
single-particle P, T-odd interaction between the proton and core ||, 

^ = ^' V ^W' < 19) 

Here it has been assumed that the proton and neutron densities are proportional to the total 
nuclear density Pa( v )] the dimensionless constant rj = jj] pp + ^r] pn . Notice that there is only 
one surviving term from the P, T-odd nucleon-nucleon interaction fll8|) ; this is because all 
other terms contain the spin of the internal nucleons for which (cr) = 0. The shape of the 
nuclear density pa and the strong potential U are known to be similar; we therefore take 
Pa( t ) — U (r). Then we can rewrite Eq. ([19]) in the following form: 

*'=^ w - ^^fr- 2 * 10 "'' 611 ' (20) 

Now it is easy to find the solution of the Schrodinger equation including the interaction W 
§: 

(H + W)i) = Ei/j , 

^ = (l + ^-V)^, (21) 

where ip is the unperturbed solution {Hip = Eip). The valence proton density is then equal 
to 

p = ^ = ,0t0 + £ V • (if) *o"ip) . (22) 
The second term gives the P, T-odd part of the density which generates a Schiff moment S 
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(23) 



where we denote r^ x as the mean-square radius of the external nucleon (in this case that of 
the proton), r 2 q is the mean-square nuclear charge radius, (cr) = tjj, and 
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i 

1+1 



I 
I 



1 + 1/2 
I - 1/2 



(24) 



where I and / are the total and orbital angular momentum of the proton, respectively. It 
should be noted that for the Schiff moment the recoil effect (the motion of the nuclear core 
around the center-of-mass) disappears due to the cancellation of its contributions to the first 
and second (screening) terms in Eq. @ ||. 

To calculate the local dipole moment, it is enough to substitute the external proton 
density ( ^2j) into the expression for the LDM (|T6| ) and perform integration using integration 
by parts. The result is 



-E 



i b x (A; + !)(* + 4) 



k + 1 
2(7+1) 



r .k+l _ 



(25) 



(Notice that for a proton in the state S1/2, the LDM is reduced to the difference of two 
approximately equal terms (r% x — r*) for all k — 2, 4, .... This makes an analytical calculation 
hopeless when trying to estimate the LDM of a nucleus with an external proton in state 
si/2, as is the case for 203 > 205 T1.) 



B. Nuclear LDM produced by core protons. Mercury moments. 

In nuclei like 199 Hg and 129 Xe the external nucleon is a neutron. It does not contribute 
to the Schiff moment directly. In Ref. [JT3] it was shown that virtual excitations of the 
core nucleons caused by a P, T-odd interaction with the external nucleon produce a Schiff 
moment which is comparable to that produced by an external proton. The actual calculation 
of the Schiff moments in |HJ was carried out numerically. In this section we perform a 
simple analytical calculation which allows us to estimate the contribution of the relativistic 
corrections ~ Z 2 a 2 to the Schiff moment. Here we follow an approach which was used in 
|i~3|l to estimate the contribution of the giant dipole resonance to the nuclear EDM. 

The expression for the local dipole moment L induced in the nuclear state |0) by the 
P, T-odd interaction W ab (0) between the nucleons a and b is 

L = (0\Mn)(n\W ab \0) + (0\W ab \n) (njLjO) (2g) 
_ (0|[g,L]|n)(nfc|0) - (0\W ab \n) (n\[H , L]|0) 

~\ (e - E n y ■ {27} 

Here H is the Hamiltonian, and [H, L] is a commutator; the LDM operator L is defined 
from Eq. ( |16D as L = (L). Now we assume that the transition strength in the sum over 
intermediate states \n) is concentrated around the excitation energy u r , and replace (E — 
E n ) 2 by uj 2 . (Note that the replacement of (E Q — E n ) by u r in Eq. ( |26| ) gives an incorrect 
result since in single-particle language there are transitions with E Q — E n = uj t and E — E n = 
—uj r .) Use of closure, J2 n \ n )( n \ = 1) gives 

L=±(0\[[H,L]M\0) . (28) 
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To calculate the commutator we assume that the motion of each nucleon in the nucleus can 
be described by the Hamiltonian H = £— + V(r). The contribution to the LDM from a 
single proton is then 

L = ^(0|(V Q L)(VJ^)|0>. (29) 

muz 

As a check of the validity of this approach for the calculation of the core contribution, 
we use this formula ( p9|) to calculate both the contribution of the external proton (which we 
can compare with Eq. (|25|) ) and the contribution of the core protons. 



1. External proton contribution 
For the external proton we can just substitute the effective potential W fll9|) into expres- 



sion (f?9j). 



G 



1 



y/2 2m 2 w r 2 



77<0|(V a L)(V a (<7- V)p A )|0>, 



(30) 



where here |0) corresponds to the state of the external proton. If we consider the nuclear 
density to be proportional to the nuclear potential U, and use the oscillator model to 
approximate the potential so that 



p A (0) 1 2 2 
H K ' U(0) 2 



then the LDM is reduced to 



2 00 b k 



_ ^ fc =i&i(^ + l)(^ + 4) LV" ' 2(1 + 1 



ti + 



k + 1 



(31) 



(32) 



Because here we are considering the LDM produced by a single proton, we can set uj r = u, 
so the factor {u 2 /uj 2 r ) = 1. Therefore, using the resonance method we can reproduce Eq. 



2. Core contribution 



In this section we use the "Schiff resonance" formalism to estimate the local dipole 
moment produced by the core protons. Again we start from Eq. (p9"D, where in this case the 
derivatives are with respect to the internal proton coordinates. Assuming that the proton 
density is proportional to the total nuclear density pa, we obtain 



G 



1 



2\f2mmuj 2 



T](0\<T- V[V a (pAV Q L)]|0) 



(33) 



where |0) is the state of the external nucleon. Because here we are considering the P,T- 



odd interaction W ap ( 15 ) between the external nucleon (proton or neutron, a) and the core 



S 



protons, the P, T-odd dimensionless constant 77 = (Z/A)r) ap . Using (|3lD to approximate the 
nuclear density, the LDM becomes 



_ _ u 2 ~ b k e£ 



r fc= 



fc 2 + 7A: + 8 / HI \ w fc + 4 fc+1 



J&i (A: + + 4) [ 2 V, 2(7 + 1 




(34) 



Before we consider the size of the relativistic corrections, let us check that this result gives 
us a reasonable value for the Schiff moment (Za — ► 0). In the approximation of a uniform 
and spherical charge distribution, r 2 = (3/5)11%, assuming that r 2 x = r 2 , and setting the 
resonance frequency for core protons u r = 2u (frequency of giant resonance), we obtain for 
the Schiff moment of 199 Hg (which has an external neutron in the state P1/2) 

S Rg « -1.6 x 10~Vpe fm 3 . (35) 



From a numerical calculation of the Schiff moment for 199 Hg performed in Ref. |TJ| the result 



S — — 1.4 x 1CT rj np e fm was obtained. This value agrees with our analytical estimate ([35]). 
Therefore it seems that the resonance method can be used for Hg to give a crude estimate 
of the size of the relativistic corrections to the Schiff moment. 

To estimate the size of the corrections, we use the approximation of a uniform and 
spherical charge distribution, r™ = (3/(n + 3))K%, and we assume that r™ x = r™, for n = 
2, 4, .... Substituting the coefficients from Appendix [A| into the expression for the LDM (p4|). 
we see that the first correction (k = 3) to the Schiff moment for Hg is about 25%. (Note 
that the second correction (k = 5) amounts to less than 10%.) Therefore the relativistic 
corrections to the Schiff moment for Hg (and we expect for other spherical nuclei) are not 
very large; for Hg we have L m 0.75S 1 . 



C. Collective LDM of an octupole deformed nucleus 

Nuclei with octupole deformation have enhanced collective Schiff moments which may 



be up to 1000 times larger than the Schiff moments of spherical nuclei HQ. In Ref. [10 
it was pointed out that the soft octupole vibration mode produces an enhancement similar 
to that of the static octupole deformation. This makes heavy atoms containing nuclei with 
collective Schiff moments attractive for future experiments searching for T-violation. 

The mechanism generating the collective Schiff moment is the following ||[7|]. In the 
"frozen" body frame the collective Schiff moment Si n t r can exist without any P, T-violation. 
However, the nucleus rotates, and this makes the expectation value of the Schiff moment 
vanish in the laboratory frame if there is no P, T-violation. (This is because the intrinsic 
Schiff moment is directed along the nuclear axis, Si ntr = Si ntr n, and in the laboratory frame 
the only possible correlation (n) oc I violates parity and time reversal.) The P, T-odd 
nuclear forces mix rotational states of opposite parity and create an average orientation of 
the nuclear axis n along the nuclear spin I, 

/ V KM , X 

\ n z) = 2q 7( /+ X ) ' ( 36 ) 

where 
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a 



E + - E_ 



(37) 



is the mixing coefficient of the opposite parity states, K — |I • n| is the absolute value of 
the projection of the nuclear spin I on the nuclear axis, M — I z , and W is the effective 
single-particle potential fll^P . The Schiff moment in the laboratory frame is 

Q - Q / \ c 2aKM focA 

3 Z — ^intT\ n z) — intr J^J _|_ ' / 

In the "frozen" body frame the surface of an axially symmetric deformed nucleus is 
described by the following expression 

R(6) = R N (l + J2PiY m (6)) . (39) 
i=i 



To keep the center-of-mass at r = we have to fix (3\ [14"|: 



ft = -3fe + . (40) 



We assume that the distributions of the protons and neutrons are the same, so the electric 
dipole moment e(r) = (since the center-of-mass of the charge distribution coincides with 
the center-of-mass) and hence there is no screening contribution to the Schiff moment. We 
also assume constant density for R < R{9). The intrinsic Schiff moment Si ntr is then [^,[7J 

Smt[ = eZ nl -L E (j+jWfta „ ez^^A , (4i) 

20" S J (21 + l)(2i + 3) "2(W35' 



where the major contribution comes from /?2/33, the product of the quadrupole (3% and oc- 
tupole /3 3 deformations. For /3 2 ~ ^3 ~ 0.1 and Z = 88 (Ra) we obtain S^tr ~ 10 efm 3 . The 
estimate of the Schiff moment in the laboratory frame gives 

eV 

S ~ aS intr ~ 0.05 e,5 2 /?3 2 ZA 2/3 r/r3- — ~ 700 x lO^e fm 3 , (42) 



where tq w 1.2 fm is the internucleon distance, E + — E- ~ 50 keV. This estimate ( [42]) is 
about 500 times larger than the Schiff moment of a spherical nucleus like Hg. Note that S 
in Eq. (fE^) is proportional to the squared octupole deformation parameter According 
to [|l(J, in nuclei with a soft octupole vibration mode ~ (0.1) 2 , i.e., this is the same as 
in nuclei with static octupole deformation. This means that a number of heavy nuclei can 
have large collective Schiff moments. 

With no screening term, it is easy to calculate the collective LDM . Use of Eq. (|TH) gives 

^4\?3(^^ 1 H i -°- 35ZV) ' (43) 

As with spherical nuclei, we see that the correction to the Schiff moment for collective nuclei 
is not very large (for Ra, Pu this correction ~ 15%). 
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V. P, T-ODD PART OF THE NUCLEAR ELECTRIC FIELD (SCHIFF FIELD) 



In this section we calculate the actual distribution of the P, T-odd component of the 
electrostatic potential y?(R) inside the nucleus (arising from the P, T-odd nucleon-nucleon 
interaction) for two models: for an external proton in a spherical nucleus and for a collective 
Schiff moment which appears in a nucleus with octupole deformation. It is found that in the 
collective case the electric field is constant and directed along the nuclear spin. This field 
distribution is also approximately correct in the spherical case when the external proton is 
in state S\/2', this is also true without the P, T-odd interaction but when the external nucleon 
(proton or neutron) possesses an intrinsic EDM and is in state Si/ 2 . 



A. P, T-odd electric field produced by a valence proton 



To calculate the P, T-odd part of the electrostatic potential tp? produced by the external 
proton, we substitute the P, T-odd perturbed external proton density (p2^) into (|6|) and 
integrate by parts, 

VPt(R) = eeV • [/ Tr^A] - ^(tr) • V J T^~^d 3 r , (44) 

where p a = ip^crip is the spin density, (cr) = tjj. Note the similarity of this expression and 
that for a P, T-odd potential produced by an external proton electric dipole moment d p (see, 
e.g., P,[12|j). In the latter case one should only replace e£ by d p (or by d n in the case of Hg 
or Xe). Note, however, that generally speaking p a = ip^cnp 7^ (a)^^ (this assumption was 
used in ||), i.e. the direction of the external nucleon EDM depends on the coordinate r. 
The separation of the spin and the coordinate variables is possible in the case of / = I + 1/2. 
Taking I = I z we obtain p a = jPm, where pu is the density of the valence proton. If we 
now assume that this density is constant within the sphere of the radius Rm, Puij) = 1 Is , 

M 

r < Rm, and Pm{t) = for r > Rm, and, similarly, the nuclear charge density p q (r) = 4 ^ , 
r < R N , and p q {r) = for r > P^v, then we obtain for the dipole term (and for R N < Rm)'- 



t 1 ~ w) R< R 



N 



<^ (1) (R) = -e£R -jUjzt-w) R n <R<R m (45) 

M R> R M . 

Thus, the P, T-odd part of the electrostatic potential is cp^(R) cx R cos 9 inside the nucleus. 
This gives us a very simple picture for the P, T-odd electric field (Schiff field): E z = — ^ oc / 
inside the nucleus. Thus, the Schiff moment gives a constant electric field along the nuclear 
spin, E oc I, and this field vanishes within the nuclear "skin" (see Fig. [I]). 

We can easily establish a relation between the P, T-odd electrostatic potential inside the 
nucleus ip^ 1 ' oc Rcos9 and the Schiff moment S. Comparing Eq. ([45]) with Eq. (|23|) (with 
tj — 1, I = 1/2) we obtain 

<p W(B) = -^^n(R-R s ) 1 (46) 

TV 
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where n(R — Rn) is a smoothed step-function Q(Rn — R), that is n(R — Rn) ~ Q(Rn — R)', 
n(R — R N ) = 1 for R < R N and ti^R—Rn) = for R > R N + 5, where 5 = Rm~Rn « Rn- 
It gives the natural generalization of the Schiff moment potential (|1|) for the case of a finite- 
size nucleus. Of course, in the general case, the radial function n(R — Rn) in the first 
harmonic of the P, T-odd potential ( |4*6| ) is more complicated (this gives some "wiggling" of 
the electric field inside the nucleus). 



B. P, T-odd electric field produced by a collective Schiff moment 

Now we wish to calculate the electrostatic potential ip^ arising due to a collective Schiff 
moment in a nucleus with octupole deformation. We use Eqs. (|39| , |40|) and assume that the 
distributions of the protons and neutrons are the same (therefore e(r) = 0, and so there is 
no screening term) and that the density for R < R(6) is constant. Calculating the integral 
in Eq. (^) for R < min R(6) gives 

yM(RH --g|L R i E C + Dftft» . (47) 



In the laboratory frame the result differs by an extra factor (n z ) fl36|). Using Eq. fl4T|) we 
can present the final result for (p^ as 



^(R) = - 1 -^^n(R-R N ) , (48) 



where n(R — Rn) = 1 for R < min R{0), and n(R — R N ) = for R > max R{6). This result 
is similar to Eq. (^6|). Thus, the collective Schiff moment produces a constant electric field 
along the nuclear spin inside the nucleus and zero field outside (Fig. p. In this case the 
width of the transition area of the nuclear surface ~ PiRn ~ 0.2R^. 

The expression ( |46| . [48|) contains the fifth power of the nuclear radius Rn which is not 



a very well-defined property in the case of finite nuclear thickness. We can avoid this 
ambiguity by using a "normalized" expression which is exact in the limit Za << 1 when 
p sp ~ R (compare with Eq. |3|): 

<^M(R) = — n(R — R N ) , B= n(R - R N )R 4 dR ^ R%/5 . (49) 

B J 



VI. THE ATOMIC EDM INDUCED IN PLUTONIUM 

In Ref. |T0] it was shown that 239 Pu has a large vibrational Schiff moment and it was 



pointed out that it is a good candidate for experiments searching for P, T-odd effects: it has 
a ground state nuclear spin 1=1/2 and it has a long half-life. The ground-state electron 
angular momentum is J = 0, and so the atomic EDM is sensitive to the nuclear Schiff 
moment; however, it corresponds to a complex electronic configuration, 5f e 7s 2 . 

In this section we perform a simple analytical estimate of the size of the atomic electric 
dipole moment induced by the nuclear Schiff moment for 239 Pu. The contribution of electrons 
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from the /-shell is small. Therefore, in a simplistic model we can consider 239 Pu to be an 
electronic analog of 199 Hg. We can then exploit the Z-dependence of the induced atomic 
EDM, as was done, e.g., in Ref. [0, to estimate the EDM of 239 Pu from the calculation of 
the atomic EDM of 199 Hg. The arguments for this simple estimate follow. 

We see from (0) that there are three factors contributing to the atomic EDM: the electric 
dipole transition amplitudes, the energy denominators, and the Schiff matrix elements. The 
first two factors are sensitive to the wavefunctions at large distances, and these are similar 
for analogous atoms. The matrix element of the Schiff moment is determined from distances 
close to the nucleus, and therefore the significant contributions come from the matrix ele- 
ments of S1/2 and p\/ 2 as well as Sy 2 and p 3 / 2 . These matrix elements strongly depend on 
the nuclear charge: they are proportional to SZ 2 Ri/ 2 and SZ 2 R 3 / 2 , respectively, where the 
relativistic enhancement factors R\/ 2 and R 3 / 2 are given by [|5| 



„ _ 47i/2 , 2ZR N \2j 1/2 -2 

1/2 ~ [r(2 7l/2 + l)] 2 ^ a B ) 

48 /2ZRn \ 71/2+73/2-3 

r(2 7 i/ 2 + l)r(27 3 /2 + 1) ^"oj 



*3 /2 = „,„., , Z,n f^y^- , (50) 



where = [(j + |) 2 — [Za) 2 \ 2 and as is the Bohr radius. Because there are twice as many 
P3/2 states as py 2 states, we use a linear combination of R\/ 2 and R 3 / 2 in the calculations, 
R sp = (R1/2 + 2i? 3 / 2 )/3. 

Our estimate for the atomic EDM induced in 239 Pu due to its nuclear Schiff moment can 
therefore be expressed in terms of the Schiff moment and atomic EDM of 199 Hg, for which 
calculations have been performed, 

4t om (Pu) = 4t om (Hg) llf 2 l SP \ Pn ■ (51) 



We use the results of Ref. |L3| for the value of the Schiff moment of Hg, Sfig = — 1.4 x 
W~ 8 r) np e fm 3 , and the atomic EDM it induces, rf at0 m(Hg) = 5.6 x I0~ 25 r] np e cm. The atomic 
structure ratio (Z 2 R sp )p u / (Z 2 R sp )i{g = 2.6. The atomic EDM for 239 Pu in terms of its Schiff 
moment is then 

4tom(Pu) = -1 x 10- 16 (S Pu /(e fm 3 ))e cm . (52) 



If we take the value for the Schiff moment of Pu from Ref. QTO | , Sp u = 400 x 10 _8 ?7 ra e fm 3 , 
then the atomic EDM is <i atom (Pu) = —4 x 10~ 22 r] n e cm, which is ~ 700 times larger than 
the atomic EDM induced in 199 Hg. 

A numerical calculation of the atomic EDMs induced in Hg, Xe, Rn, Ra, and Pu is 
underway. 
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APPENDIX A: ELECTRON TRANSITION DENSITY FOR S-P 1/2 AND S-P 3/2 

To calculate the electron wavef unctions inside the nucleus R < Rn we assume that the 
nuclear charge is uniformly distributed about a sphere. This charge distribution corresponds 
to the harmonic-oscillator potential 



V 



Za /3 



R 



N 



-X 



(Al) 



where we have set x = R/Rn- By solving the radial Dirac equations for an electron in states 
s i/2, P1/2, and p 3 /2 moving in this potential we obtain the radial wavefunctions (see Eq. 
for definition) for s: 



f 8 (R) = f s (0)\l - -Z 2 a 2 x 2 (l - —x 2 + —x 4 - —Z 2 a 2 x 2 + 
Js\ 1 Jsk n g V 15 45 80 

1 / 1 9 

— R^mZax 2 \\ x 2 Z 2 a 2 x 2 + ... 

9. V m 40 



g s (R) = f s (Q) 



2 v 10 40 

Zaxil x 2 Z 2 a 2 x 2 + ... 

2 V 5 40 



+— R N mZ 2 a 2 x 3 (l - —x 2 + —x 4 - —Z 2 o?x 2 + .. 
20 V 42 54 56 



for p 1/2 : 



fpyA R ) = 9 Pl/2 (0)[^Zax(l - ^x 2 - ^-Z 2 a 2 x 2 + .. 



+ -R N mx\l - — Z 2 a 2 x 2 (l - —x 2 + 



324 



27 
224 



Z 2 a 2 x 2 + ...)} 



9*,M = g^Ml 1 ~ \zw{i - ±x 2 + lx 4 - JLzw + ...) 



— RfjmZax 2 {\ x 2 — ^-Z 2 a 2 x 2 

2 V 6 40 



and for p 3/2 : 



fpz/ 2 ( R ) = {fp 3/2 /x) R=0 x[l - — Z 2 a 2 x 2 (l - -x 2 + - 



() "* * " 1 " 5 - x * - JLz 2 a 2 x 2 + ... 



59 



112 



9 P3/2 ( R ) = -Za{f P3/2 /x) R=0 x 



A _ l x 2 _ IL Z 2 a 2 x 2 + .. 
10 14 560 



Here / s (0), g Pl/2 (0), and (f Ps/2 /x)R=o are the s, py 2 , and (p 3 / 2 /x) radial wavefunctions at 
zero, and m is the electron mass. The terms included into the radial wavefunctions above 
are such that the radial transition densities U sp = f s f p + g s g p for s-py 2 and s-p 3 / 2 include 
all corrections of order Z 2 a 2 and the lowest correction of order Z 4 a , 

U sPl/2 = \fs^)g Pl/ MmR N x{l - ?ZW(l - ^-x 2 + ^-x 4 ) + ^~Z 4 a 4 x 4 + ...} (A2) 



SP3/2 



W)(f P3/2 /x) R=Q x{: 



20' 



-Z 2 a 2 x 2 (l 



14 
69 

315' 



:X 



560 

—x 4 ) + —Z 4 a 4 x 4 + ...\ . (A3) 
63 / 2800 J v ' 
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It is seen by direct substitution of the transition densities (|A2| , |A3|) into the LDM expres- 
sions (^,^4],[E|) that it is sufficient to include in U sp just the first correction. These terms 
(oc Z 2 a 2 /R 2 N ) correspond to the coefficient 63 in expression (0). They give corrections 
to the Schiff moments of ~ 15 — 25%. The remaining terms in (|A2|JA3|) correct the Schiff 
moments by a few percent. 
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FIGURES 




FIG. 1. Constant electric field E inside the nucleus produced by the Schiff moment. E is 
directed along the nuclear spin I. 



